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The time evolution of temperatures of anisotropic nanoparticles in two and three-body systems 
are simulated for various relative orientations. Nanoparticles are immersed in a thermal bath at 
constant temperature. It is shown that in two-body systems, the relative orientation of nanoparti¬ 
cles could drastically affect the dynamics of temperature evolution and thermalization time scale. 
Moreover, in some configurations, the temperature difference in initial state has a minor effect on 
the dynamics of temperatures. In three-body systems, the orientation of the third nanoparticle 
influences the temperature dynamics, which allows one to control the thermalization time scales 
between anisotropic nanoparticles. Also, in addition to previously known contribution of the small¬ 
est distance between isotropic nanoparticles on the thermalization time scales, it is shown that the 
nanoparticles’ orientations are more important in some particular arrangements. 
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I. INTRODUCTION 

It is well-known that the radiative heat transfer in¬ 
tensity can be sufficiently enhanced when the separation 
between objects is much smaller than the characteris¬ 
tic wavelength of thermal emission. This phenomenon is 
demonstrated experimentally for plane-plane and sphere- 
plane geometry and it is shown that the heat flow at 
nanoscale distances is greater than the black body radi¬ 
ation law posed by Stefan Boltzmanic— . Several theo¬ 
ries have been proposed to explain the heat transfer be¬ 
tween simple particle geometries at nanoscale distances 
and have been widely used for two-body systems in¬ 
cluding parallel slabs^—, a particle in front of a pla¬ 
nar surface^ ' 10 : 11 ; two spheres^— , and two anisotropic 
objects^— . There are series of theoretical studies in 
recent years concerning the many body effects in the ra¬ 
diative heat transfer— 5 . It is shown in many publi¬ 
cations that the heat flux exchange between two objects 
can be remarkably enhanced in a three-body system with 
respect to a two-body system, thanks to many-body in¬ 
teractions. In [25j , Ban-Abdallah et al. used the fluc- 
tuational electrodynamics theorem to develop a theoret¬ 
ical formalism for radiative heat transfer between spher¬ 
ical nanoparticles in a many-body system in the dipole 
limit. They have shown a configurational resonance in 
heat exchange between nanoparticles due to many-body 
interactions. Similarly, the exaltation of the heat trans¬ 
fer in a three-body system of parallel slab is reported 
in (26(. There are very few studies in the heat trans¬ 
fer between anisotropic materials. The effect of shape, 
orientation and polarizability tensor in the heat trans¬ 
fer between anisotropic objects allows for a large free¬ 
dom of tunability. Early research on the heat exchange 
between anisotropic objects demonstrate that the rela¬ 
tive orientation and polarizability tensors offer additional 
mechanism to tune the heat transfer, which is of inter¬ 
est for the design on new materials capable of tailor¬ 
ing heat propagation and applications in thermal sen¬ 
sors. Huth et al. studied the shape-dependence of heat 


transfer between a shperoidal nanoparticle and a pla¬ 
nar surfaces for various distances— Compared to the 
sphere-plane geometry, they have shown that it is pos¬ 
sible to control the heat transport in the system by the 
shape. Biehs et al. showed that the heat exchange be¬ 
tween two anisotropic media can be significantly larger 
than those traditionally measured between two isotropic 
materials— The strong dependence of heat exchange be¬ 
tween two anisotropic nanoparticles on their shapes and 
relative orientations has been recently reported and inter¬ 
preted as a heat transfer switch— In our recent work—, 
we used the fluctuational electrodynamics to introduce a 
formalism to take into account the change in particle’s 
shape and more generally the polarizability tensors in a 
many-body system. We showed the heat exchange in a 
system of anisotropic nanoparticles can be tuned drasti¬ 
cally depending on relative orientations and polarizabil¬ 
ity tensors. 

Most of the theoretical efforts in the investigation of 
radiative heat transfer have been restricted to the steady 
state regime. The issue of radiative thermalization and 
dynamics of temperatures has been addressed in few pa¬ 
pers and very recently some concrete ideas have been put 
forward. Tschkin et al. have investigated the radiative 
cooling of a spherical nanoparticle close to an interface 
with fixed temperature— They have shown that the 
thermal relaxation time of the nanoparticle varies rapidly 
for distances smaller than the thermal wavelength and 
is very sensitive to the initial temperature differences. 
Dyakov, et al. studied the propagation of thermal radi¬ 
ation in a two layered structures. In [29], they assumed 
that the temperature of one of the plates is changed due 
to the radiative heat transfer from another plate at con¬ 
stant temperature which acts as a heat source. They have 
shown that the thermalization time scale between the 
plates strongly depends on the separation distance. The 
dynamics of heat transfer between isotropic nanoparticles 
is well treated by Messina and his co-workers with con¬ 
sideration of many body effects— Using a purely fluctu¬ 
ational electrodynamics approach, they studied the time 
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evolution of temperatures of spherical nanoparticles with 
respect to several geometrical arrangements. It was sug¬ 
gested that many-body interactions might lead to rapid 
thermalization. They concluded that the smallest dis¬ 
tance between nanoparticles in a system determines the 
time scale on which the heat flux is exchanged between 
the particles. 

In this paper, we make a step forward in simulating the 
dynamics of heat transfer between anisotropic nanopar¬ 
ticles which are immersed inside a thermal bath. To do 
this, we use the fluctuational electrodynamics theorem 
which relates the power spectral density of fluctuating 
charge density to the local temperature and frequency 
dependent relative dielectric response of nanoparticle^. 
Using this theorem and properties of the dyadic Green’s 
function as well as polarizability tensors, allows study¬ 
ing the phenomenon of heat transfer between anisotropic 
nanoparticles at different temperatures. For the sake of 
simplicity, we will restrict ourselves to nanoparticles with 
spheroidal shape only which made of silicon carbid (SiC). 
We are interested in the effects of anisotropy, i.e., the 
shape and the orientation of shperoidal nanoparticles, on 
the dynamics of temperatures as well as the thermal¬ 
ization time scales in two and three-body systems. We 
explain our findings by studying the heat transfer as a 
function of relative orientation of nanoparticles. 

The structure of the paper is as follows. Our study 
of dynamics of temperatures in system of anisotropic 
nanoparticles begins in section. ITT] But before digging in, 
we will discuss briefly the basic ideas of heat exchange 
between anisotropic nanoparticles in presence of an ex¬ 
ternal thermal bath. In Sec. mu we consider dynamics 
of heat transfer between two spheroidal nanoparticles of 
silicon carbid (SiC) immersed in a thermal bath. We 
studied the dependence of the temperatures’ dynamics 
on the relative orientation of nanoparticles. To give an 
intuitive realization, the phase portraits of the system 
are also presented for various configurations. Section HVl 
contains the results of temperature dynamics in three- 
body system of anisotropic nanoparticles. This section is 
divided into four subsections which the dependences of 
thermalization time scales on the orientation of the third 
particle are discussed in Secs. IIV A) to IIV Dl Finally, we 
draw our conclusion in Sec. m 


II. DYNAMICS STUDY OF HEAT EXCHANGE 
BETWEEN ANISOTROPIC PARTICLES 

Before we go through the dynamics of heat transfer 
between anisotropic particles, let us focus on the prob¬ 
lem of radiative heat exchange in a system of anisotropic 
nanoparticles without considering the temporal evolution 
of temperatures. Our theoretical framework for heat ex¬ 
change between anisotropic nanoparticles has been de¬ 
scribed in previous publication 27] and is briefly summa¬ 
rized here. We considered an energy exchange problem 
between N anisotropic nanoparticles identified for con¬ 


venience by indexes i = 1, 2, • ■ • ,N with positions (rj), 
polarizability tensors (eu), and temperatures (Tj) . In or¬ 
der to take into account the radiation impinging on such a 
system from the outside, we assumed that this collection 
is immersed in a thermal bath at a constant tempera¬ 
ture (T;,). Based on the assumption that nanoparticles 
behave as fluctuating dipoles at different temperatures, 
the fluctuation-dissipation theorem provides the energy 
which is dissipated in each nanoparticles. The nanoparti¬ 
cles themselves interact through the radiation fields and 
on the other side with the thermal bath and the total ra¬ 
diation power dissipated in each particle were obtained 
as: 
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Where J)(Tj) denotes the so-called self-emission by i-th 
particle corresponding to the heat lost by its radiation 
in presence of other particles in the system, ^Fij^iiTj) 
is the radiative heating of the i-th particle due to radi¬ 
ation of the j-th one, and finally JF\ 3] t (7b) is a portion 
of power which is absorbed in the *-th particle due to 
the radiation of the thermal bath. In these equations, 
0(w,T) = huj[ 1 + 2n(w,T)]/2 and n(ui,T) is the Bose- 
Einstein energy distribution function and we have intro¬ 
duced x = a + fc 2 <iGjad. Here G,;, = G ji are the elec¬ 
tric Dyadic Green functions, Go = and K = oj/c. 

The physical properties of nanoparticles and the mathe¬ 
matical effects of many-body interactions are contained 
in the coupling matrixes A, B, C and D. For the special 
case of three anisotropic nanoparticles, we find 
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The frequency-dependent relative dielectric permittivity 
of SiC is well described by lorentz rnodel^: 


e(w) 



9 2 

UJ L — 

— uj 2 — iTto 


( 4 ) 


with uj l = 969 cm -1 , u>t = 793 cm -1 , T = 4.76 cm -1 
and £oo = 6.7. This relative permittivity is used to cal¬ 
culate the dressed polarizability tensors of nanoparticles. 
In the case of small isotropic non-spherical particles (or 
different extensions, including anisotropic spherical par¬ 
ticles or anisotropic non-spherical particles), the dressed 
polarizability tensor can be expressed as^: 


r p i -i 
6l% = cx-hl 1 - i—OLh > , 


( 5 ) 


with 


ct h = i>i(i - i)< i + (L - fc 2 M)(£ - i) 


( 6 ) 


where L and M are the singular and non-singular part 
of the depolarization dyadio^i. For special case of small 
prolate spheroidal nanoparticles, the polarizability tensor 
is diagonal in the principal-axis system of each particle 
and the components of the polarizability tensor along and 
across the rotation axis of a spheroidal nanoparticle are 
given by: 


4 2 £ ( UJ ) ~ 1 

a h,P — . r , , , 

3 1 + Lp[e(w) - 1J 


,0 8 = b,s), (7) 


where Rb and R s are the major and minor axeses of 
spheroid and the depolarization factors take the form^: 
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Now that we have established the formalism of heat ex¬ 
change between anisotropic nanoparticles, we can present 
a framework for dynamics of temperatures. The time 
evolution of the temperatures in system is a solution of 
the nonlinear coupled system of differential equation: 

d T. 

Pi c i v i ~oT” = 'P i (Ti,T 2 ,--- , T N ,T b ) (z = 1,2,--- , IV). 
at 

. (10) 

where Vi is the volume of the particle with mass density pi 
and heat capacity Cj. The temperature-dependent heat 
capacity of SiC is calculated as follows^: 

c = 1.267+0.049 x 10 _3 T—1.227 xl0 5 T _2 +0.205 x 10 8 T -3 

( 11 ) 


We note that the right-hand side of Eq. m has N non¬ 
linear terms with a strong nonlinearity in temperatures. 
We need the initial values of temperatures as well as the 
bath’s temperature to determine the solution uniquely. 
The solutions of this equation could be visualized as tra¬ 
jectories flowing through an N-dimensional phase space 
with coordinates (Tj,-- - , Tjv). We will discuss the dy¬ 
namics of temperatures, and will show that the qualita¬ 
tive structure of the temperatures’ dynamics can change 
as the relative orientation of particles is varied. 


III. DYNAMICS OF HEAT TRANSFER IN 
TWO-BODY SYSTEM 


In this section we study the temperature evolution in 
two-body systems and visualize the dependence of tem¬ 
perature dynamics on the nanoparticles relative orien¬ 
tation. The system consists of two prolate spheroidal 
nanoparticles (major to minor axis ratio, r) = 4) with a 
fixed center to center separation, d = 800 nm. The parti¬ 
cles are heated up to initial temperatures Tj(0) and T 2 (0) 
and immersed in a thermal bath at fixed temperature Tb- 
The nanoparticles themselves interact with each other 
and also with the thermal bath. Our aim is to calculate 
the particles’ temperatures, Tj and Tj, as a function of 
time. So, Eq. (1) is used to calculate the total power dis¬ 
sipated in each particle and the dynamics of temperatures 
are calculated according to the energy balance equation, 
Eq. (1101) . Since, the thermalization rate in each parti¬ 
cle due to thermal conduction by heat carriers is much 
larger than the heat transfer rate between the particle 
and the rest of the system (including other particles and 
thermal bath) via radiation, it ensures that the energy 
fluxes supported by heat carrier inside each particle is 
large enough to treat the temperatures as homogeneous 
inside each particle during the dynamical evolution of the 
system. 

In order to show the relative orientation effect on the 
dynamics of heat transfer, particles are arranged in the 
system with one of the main arrangements, namely the 
tip-tip, side-side, tip-side (side-tip), and cross configura¬ 
tions. Figure Q shows solution to Eq. m for these 
types of configurations. The initial temperatures are 
Ti = 350 K and T 2 = 300 K and particles are sepa¬ 
rated d = 800 nm apart in a thermal bath at Tb = 300 
K. The fixed point of Eq. (flOl) corresponds to the ther¬ 
mal equilibrium state of the system. In the presence of 
thermal bath (at constant temperature) and absence of 
any other external sources, the particles would thermal- 
ize to the temperature of the bath. On the other hand, 
there is a fixed point at bath’s temperature, regardless of 
the nanoparticles’ characteristics or initial temperatures. 
We would expect no change in temperatures in thermal 
equilibrium and from Fig. © it is clear that it takes ap¬ 
proximately 0.1 second for particles to thermalize with 
bath ( i.e., Ti j2 —>• 300 K) which is independent from the 
relative orientation of the nanoparticles. However, it can 
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FIG. 1. (Color online) Time evolution of the temperatures in 
a two-body system for (a) tip-tip, (b) side-side, (c) cross and 
(d) tip-side (side-tip) configurations. The distance between 
particles is d = 800 nm and initial temperatures are Ti = 350 
K, T 2 = 300 K and thermal bath is at fixed temperature 
Tb = 300 K. Nanoparticles are prolate spheroid with major 
to minor aspect ratio g = Rb/R s = 4. The red (solid) line 
corresponds to the hotter particle while the blue (dashed) line 
is for cold one. 


be seen that there exist two types of time scales in the 
curves; one corresponds to the thermalization between 
two particles (smaller time scale) and the other for ther¬ 
malization of particles with the bath (larger time scale). 
The former is mainly due to the intra-particle heat ex¬ 
change while the later is caused by the radiative heating 
of nanoparticles by the thermal bath. 

In the case of tip-tip configuration, the first thermal¬ 
ization is ~ 1 order of magnitude faster than the second 
thermalization. For the side-side configuration, the early 
thermalization time scale increased slightly and the first 
thermalization temperature tends to T\ j,. This is because 
the coupling between particles in this configuration is 
small in comparison with the tip-tip case and the heat 
transported from the particle at larger temperature to 
the colder one could slightly compensate the decrease of 
the temperature of the colder one due to its thermal emis¬ 
sion. With decrease of ordering in the system, by going 
from tip-tip or side-side to cross or tip-side configura¬ 
tions, the heat exchange between particles drops several 
orders of magnitude^. As a result, the rate of energy 
lost in the cooler particle is more than received from the 
hotter one and it seems that particles are thermalize by 
the bath independently. 

The phase portraits for different configurations of par¬ 
ticles in the two-body system are shown in Fig. ©. Each 
of these configurations has a stable fixed point at T* = Tt, 
(Solid blue dot) and all trajectories approach it as t —> 00 . 
However, the direction of approach depends on the initial 
temperature difference as well as the relative orientation 
of the particles. Moreover, the thickness of the trajecto- 
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FIG. 2. Phase portrait of the two-body system immersed 
in a thermal bath at temperature Tb=300 K for (a) tip-tip, 
(b) side-side, (c) cross and (d) tip-side configurations. The 
thickness of the trajectories show the thermalization rate. 

ries shows the thermalization rate which is faster for large 
temperature differences and the phase point slows down 
as reaches the proximity of the thermalization point. 

In the case of tip-tip configuration, as shown in 
Fig. (J2Jt), all trajectories slam straight onto the slower 
direction (here, the Ti = T 2 direction) which corresponds 
to intra-particle thermalization regime, then slowly ooze 
along this curve until they reach a fixed point. The point 
at which the trajectory reaches this line is the intra¬ 
particles’ thermalization temperature. The approach to 
the fixed point is tangent to the Tj = T 2 direction which 
the thermalization by bath occurs. On the other hand, 
if we look backwards along a trajectory (t —> — 00 ), then 
the trajectories all become parallel to the faster decaying 
direction (here, the T\ + T 2 — 2T, = 0 direction). 

As the relative orientation changes to side-side con¬ 
figuration, the coupling between the nanoparticles de¬ 
creases an order of magnitude in comparison to tip-tip 
configuration^. As a result, the time scale for thermal¬ 
ization of the particles with each other increases and the 
decay to the slower direction decreased. The trajectories 
shows faster decay along the T) = T 2 direction than pre¬ 
vious configuration and approaches the fixed point mono- 
tonically as t 00 , [see Fig. (j2j>)]. Furthermore, similar 
to the tip-tip configuration, intra-particle thermalization 
temperatures depend strongly on the initial temperatures 
of the particles. 

In the case of cross and tip-side configurations, 
Figs. m and ©1) respectively, the coupling between 
nanoparticles is very week in comparison to the previ¬ 
ous configurations. In these cases, the thermalization by 
thermal bath is more effective and T* = T, is a sym- 
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metrical node or star shaped fixed point of the system 
which represents an attractor that the flows are toward it. 
Since the decay rates of temperatures are approximately 
the same, all trajectories are approximately straight lines 
through the fixed point. Moreover, the time it takes for 
particles 1 and 2 to thermalize with each other does not 
appreciably depends on their initial temperatures and is 
of order of thermalization with bath. From these phase 
portraits it becomes apparent that both the temperature 
differences and relative orientation of the particles deter¬ 
mine the time scale on which the heat flux is exchanged 
between anisotropic nanoparticles. 


IV. DYNAMICS OF HEAT TRANSFER IN 
THREE-BODY SYSTEMS 


Let us now calculate the time evolution of the tem¬ 
peratures in a three-body system of anisotropic nanopar¬ 
ticles. For the sake of simplicity, consider a two-body 
system at which particles 1 and 2 are oriented in one of 
the configuration discussed in previous section. Now, a 
third nanoparticle labeled with i = 3 with the same phys¬ 
ical properties, is placed between nanoparticles 1 and 
2 [see for example Fig. dUr)]. In such a system, each 
of the nanoparticles exchanges heat with the two others 
and the bath. Once again the dynamics of temperatures 
obey Eq. (flUl) and we will show that the system under¬ 
goes different thermal dynamics depending on the orien¬ 
tation of the third particle. Intuitively we know that the 
coupling between 1 and 2 affects by the presence of the 
third particle and leads to an enhancement (reduction) 
of the radiative heat exchange, therefore, to shortening 
(enlargement) of time scale of intra-particle thermaliza¬ 
tion. Moreover, addition of a new particle would alter 
the coupling between each particle and the thermal bath 
and apart from the intra-particle interaction it would af¬ 
fect the thermalization time scales. From a dynamical 
point of view, addition of the third particle increases the 
phase space dimension which might allows for more pos¬ 
sible dynamical properties and we may encounter some 
new phenomena. As an example, adding a particle in¬ 
troduces new intra-particle thermalization time scales to 
the problem which depends on initial temperatures. 

To compare the results with the previous section, we 
assume that particles 1 and 2 are arranged in either tip- 
tip, side-side, cross or tip-side configuration and the evo¬ 
lution of the temperatures are discussed for a particular 
orientation of the third particle in each of these types. 
The time evolution of the temperatures are calculated for 
the initial temperatures Tj = 350 K and X 2 = I 3 = 300 
K, and the thermal bath is maintained at fixed temper¬ 
ature Tf, = 300 K. 



FIG. 3. (Color online) Dynamics of temperatures in a three- 
body system. Particles 1 and 2 are oriented in a tip-tip config¬ 
uration and a third particle is placed between them to form a 
(a) tip-tip-tip and (b) tip-side-tip configuration. The separa¬ 
tion between 1 and 2 is d = 800 nm as the two body case and 
initial temperatures are Ti = 350 K, T 2 = 300 K and T 3 = 300 
K, respectively. The collection is immersed in a thermal bath 
at Tb ~ 300 K. The red (solid) line corresponds to the hotter 
particle (the first particle), the blue (dashed) line is for cold 
one (the second particle) and the green (long-dashed) line is 
for the third particle. 


A. Three-body system: particles 1 and 2 in tip-tip 
configuration 

Let us first consider the thermal evolution problem in 
a system where particles 1 and 2 form a tip-tip configura¬ 
tion. Figures [3Ji and[ 3 j 3 represents the time evolution of 
the temperatures calculated for two particular arrange¬ 
ments, namely, tip-tip-tip and tip-side-tip, respectively. 

In Fig. J3]i), nanoparticles are aligned along the longer 
spheroid axis and the intra-particle couplings are very 
strong. Surprisingly, it can be seen that the presence 
of the third particle reduces the time scale for thermal¬ 
ization between 1 and 2 of about an order of magni¬ 
tude in comparison with the two-body system. More¬ 
over, the heat exchange between 1 and 3 is enough to 
let the third particle thermalize by first one even slightly 
faster than the time scale of thermalization of 1 and 2 . 
This is explained by the fact that the coupling between 
1 and 3 is larger than the coupling between 1 and 2 
due to separation distances. After such an intra-particle 
thermalization, it takes a long time for the phase point 
(7i = T 2 = T 3 ) to pass through such a bottleneck state be¬ 
fore making a transition to the final thermalization state 
of the system. Since the coupling with the thermal bath 
does not change significantly by the presence of the third 
particle, the time it takes for the system to reach the 
steady state is as the same order in the two-body sys¬ 
tem. 

As the third particle rotates to form a tip-side-tip con- 










6 





Time (s) 

FIG. 4. (Color online) Dynamics of temperatures in a three- 
body system. Particles 1 and 2 are oriented in a side-side 
configuration and a third particle placed between them to 
form a (a) side-side-side, (b) side-tip-side, and (c) cross-cross 
configuration. 

figuration, Fig. ©), the coupling between 1 and 2 drops 
an order of magnitude in comparison with the previous 
case^l. The tip-side coupling between 3 and 1 (or 2) is 
not large enough and it has a trivial role in the dynamics 
of temperatures of particle 1 and 2. As a result, the tem¬ 
perature evolutions of particles 1 and 2 are similar to the 
two-body case [see Fig. ©>)]■ Moreover, particle 3 is not 
affected by particle 1, as much as particle 2, even though 
it is closer to it and has the same initial temperature 
as particle 2. Basically, this corresponds to the nature 
of couplings between anisotropic nanoparticles. On the 
other hand, while the smallest distances between parti¬ 
cles determines the thermalization time scales in systems 
with isotropic particles^, the orientations (or more gen¬ 
erally, the polarizability tensors) play a crucial role when 
we are dealing with anisotropic particles. 


B. Three-body system: particles 1 and 2 in 
side-side configuration 

There exist three particular orientations in adding a 
third particle to a two-body system with side-side con¬ 
figuration. Figure ©,) shows the temperature evolution 
for side-side-side configuration. Since the relative orien¬ 
tation of each couple of particles in the system are iden¬ 
tical (side-side form), the intra-particle distances deter- 




Time (s) 

FIG. 5. (Color online) Dynamics of temperatures in a three- 
body system. Particles 1 and 2 are oriented in a cross like 
configuration and a third particle placed between them to 
form a (a) cross-side-side, (b) side-tip-side, and (c) side-side- 
cross configuration. 

mines the time scales on which the heat flux is exchanged 
between the particles. In this case, particle 3 is closer to 
particle 1 and as a result it is heated faster than particle 
2. Once again, the shift in the timescales can be at¬ 
tributed to the presence of the third particle which acts 
as a bridge for indirect heat transfer between particles 1 
to 2. 

The results of the temperature evolution of particles 
for other orientation of the third particle are shown in 
Figs. ©:>) and ©), for side-tip-side and cross-cross con¬ 
figurations, respectively. The orientation of the third par¬ 
ticle does not affect the intra-particle coupling between 
particles 1 and 2 significantly. Moreover, particle 3 is 
heated slightly more by particle 1 in the cross-cross con¬ 
figuration before the final thermalization with the ther¬ 
mal bath. 


C. Three-body system: particles 1 and 2 in cross 
configuration 

The results for dynamical evolution of temperatures in 
a three-body system in which particles 1 and 2 oriented 
in cross like form are illustrated in Figs. ©i)-©). As 
can be seen in Fig. ©t), the stronger coupling between 
particle 3 and 2 results in a similar thermal evolution 
of these particles toward the steady state of the system. 
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FIG. 6. (Color online) Dynamics of temperatures in a three- 
body system. Particles 1 and 2 are oriented in a tip-side 
configuration and a third particle placed between them to 
form a (a) tip-side-side, (b) tip-tip-side, (c) tip-side-cross, and 
(d) tip-oblique-side configuration. 

In figure. [5 }d, the tip-side coupling of particle 3 with 1 
(and 2) leaves the dynamics of temperatures unchanged. 
On the contrary, as we can see in Fig. the coupling 
between particles 3 and 1 is more effective and results 
in a faster thermalization between them. This coupling 
causes a drastic change in the thermal evolution of the 
first particle. 

D. Three-body system: particles 1 and 2 in 
tip-side configuration 

Now, consider a system in which particles 1 and 2 are 
arranged in a tip-side configuration. By adding a third 
particle, with one of the three mentioned orientations, 
three possible configurations are formed [see Figs. (jGja,)- 
Efc)] . It can be seen in these figures that the temperature 
evolution of particle 2 is not influenced by the presence of 
particle 3. However the intra-particle thermalization be¬ 
tween particles 1 and 3 is very sensitive to the orientation 


of the third particle and specially pronounced in a tip- 
tip-side configuration shown in Fig. d^Jo) . The more in¬ 
teresting case is for configuration in which the third par¬ 
ticle rotates by angle 45 degree in the plane of particles 
1 and 2, namely tip-oblique-side configuration as shown 
in Fig. ©1). Upon inspection of low coupling between 
particles 1 and 2 observed in previous configurations, we 
find that the temperature evolution in tip-oblique-side 
is largely decided by three-body interactions. Here, the 
third particle appears to remove energy from the first 
particle and re-inject almost all of it to the second one. 
One can see that the temperature evolutions are drasti¬ 
cally affected by the orientation of the third particle and 
the thermalization between particles happens very fast. 

V. CONCLUSION 

Using the method of fluctuational electrodynamics in 
conjunction with the properties of the dyadic Green’s 
function as well as polarizability tensors, we have inves¬ 
tigated the dynamics of temperatures in two (and three) 
body systems of anisotropic nanoparticles. Nanoparticles 
are immersed in a thermal bath at constant temperature 
and the dynamics of temperatures are simulated with re¬ 
spect to the relative orientation of nanoparticles. In the 
case of two-body system we showed the drastic effect of 
the relative orientation of nanoparticles on the dynam¬ 
ics of temperatures in the phase space of the system. In 
the case of three-body systems, we have found that the 
orientation of the third nanoparticle influences the dy¬ 
namics of temperatures in comparison to the two-body 
systems. Also, in addition to the aforementioned works 
on the sensitivity of the heat transfer on the separation 
distances, it turns out that the dependence of the heat 
transfer on the shape and relative orientation of particles 
can be much larger than the smallest distances in some 
cases. Although we have considered the heat transfer 
dynamics for spheroidal particles, our study reveals the 
hidden characteristics of heat transfer between arbitrary 
shaped (and/or polarizability tensor) objects. We believe 
that our findings are essential for understanding the heat 
transfer in anisotropic systems and the physical results 
of this work relevant to control the thermalization time 
scales in nanoscale devices. 
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